Amorphous nucleation precursor in highly nonequihbrium fluids 
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Dynamical density functional simulations reveal structural aspects of crystal nucleation in under- 
cooled liquids: the first appearing solid is amorphous, which promotes the nucleation of bcc crystals, 
but suppresses the appearance of the fee and hep phases. These findings are associated with features 
of the effective interaction potential deduced from the amorphous structure. 
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Mounting evidence indicates that the classical picture 
of crystal nucleation, which considers "heterophase" fluc- 
tuations of only the stable phase, is oversimplified. Early 
analysis by Alexander and McTague suggests preference 
for bcc freezing in simple liquids [1 . Atomistic simu- 
lations for the Lennard- Jones (LJ) system have verified 
that small heterophase fluctuations have the metastable 
(MS) bcc structure, and even larger clusters of the sta- 
ble (S) fee structure have a bcc interface layer [2 , while 
the ratio of the two phases can be tuned by changing 
the pressure [3 . Composite bcc/fcc nuclei have also been 
predicted by continuum models [4 . Two-stage nucleation 
has been reported in systems that have a metastable criti- 
cal point in the undercooled liquid (including solutions of 
globular proteins [5 and eutectic alloys 0); the appear- 
ance of the crystalline phase is assisted by liquid droplets, 
whose formation precedes/helps crystal nucleation [7 . 
Recent studies indicate a similar behavior in simple liq- 
uids such as the LJ [8 or hard-sphere (HS) [9^ fluids, 
where dense liquid/amorphous precursor assists crystal 
nucleation. Analogous behavior has been reported for 
colloidal systems in 2D [10 and 3D pJJ. Brownian Dy- 
namics studies for the HS system p!2l show the evolu- 
tion of medium range crystalline order during the pre- 
nucleation stages. Liquid mediated crystal- amorphous 
and crystal-crystal transitions have also been predicted 
[13] . These findings imply that the nucleation precursors 
are fairly common. A deeper understanding of nucle- 
ation pathways requires a systematic study of a system, 
in which amorphous and crystalline structures compete 
during solidification. 

Such a system is defined by the single-mode phase-field 
crystal (IM-PFC) model of Elder et at. [14 , a simple 
dynamical density functional theory, which has bcc, fee, 
and hep stability domains [15], and the appearance of 
amorphous phase and two-step nucleation has also been 
reported [16 . A two-mode extension of the model by 
Wu et al. (2M-PFC) has been designed to promote fee 
solidification [17 , whereas with a specific choice of model 



parameters the IM-PFC model can also be recovered. 

Herein, we address crystal nucleation in PEG models 
interpolating between the IM-PEC and 2M-PEG limits. 

First, we re-cast the free energy of the PFG models 
in terms of A = Ri/{1 + Ri) G [0,1] {Ri is the rela- 
tive strength of the first- and second-mode contributions 
[17 ), a parameter that can be used to interpolate be- 
tween the 2M-PFC (A = 0) and IM-PFG (A = 1) models: 

:F= /dr{|[e + (l + VY 

4 (1) 

x(A + {l-A}{Q? + vn^)]V^+^}, 

where ij) oc {p — P^l"^)/ p]^"^ is the scaled density difference 
relative to the reference liquid of particle density p^^*^. 
The reduced temperature e can be related to the bulk 
moduli of the fluid and the crystal, whereas Qi = qi/qo 
(= 2/3^/2 for fee [TT) is the ratio of the wave numbers 
corresponding to the two modes. 

The respective dimensionless Euler-Lagrange equation 
(ELE) and equation of motion (EDM) read as |^ = 

(f . and H = V^ff + C, respectively, where ff 

denotes the functional derivative of with respect to 
t/j^ and r is the dimensionless time. The RHS of the 
ELE is taken at the far-field value i/jq (homogeneous liq- 
uid). In the EDM, the fluctuations are represented by a 
colored Gaussian noise ( of correlator (C(r, r)C(r', r')) = 
— aV^^(|r — r'l, <j)(5(r — r'), where a is the noise strength, 
and g{\r — r'l^cr) a high frequency cutoff function [18] 
for wavelengths shorter than the interatomic spacing (a) . 
Due to the overdamped conservative dynamics the EDM 
realizes, the PFG models defined so are suitable for de- 
scribing crystalline colloidal aggregation [191 EO] • These 
equations have been solved numerically [21 on rectan- 
gular grids of typical size of 512 x 256 x 256 (ELE) and 
256 X 256 X 256 (EDM), assuming a periodic boundary 
condition. The ELE has been used to determine the 
phase diagram, the driving force, the nucleation barrier. 
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FIG. 1: (color online) Two-step nucleation in the IM-PFC 
model at e = —0.1667 and ipo = —0.25. Left: Snap- 
shots of the density distribution taken at the dimensionless 
times T = 57.74t. Spheres of diameter of the interparticle 
distance centered on density peaks higher than a threshold 
(0.15) are shown that are colored red if G [0.02,0.07] and 
qe G [0.48,0.52] (bcc-like), and white otherwise. Right: pop- 
ulation distribution of qe (histogram painted similarly), and 
the time-dependence of the fraction X of bcc-like neighbor- 
hoods (solid line). Note the nucleation of amorphous clusters 
and the formation of amorphous grain boundaries [23^. 



and the coexistence properties, including the densities 
and the solid-liquid interface free energy (as described in 
[15 ), while the EOM has been applied to simulate nu- 
cleation. Owing to the effect of noise on the free energy, 
the results from the two approaches converge for ^ ^ 0. 

The results of the nucleation studies performed solv- 
ing the EOM under condition described in [22] are sum- 
marized in Figs, [l] and [2] In the case of the IM-PFC 
model, in a large part of the bcc stability domain we 
have observed two-step nucleation starting with forma- 
tion of amorphous clusters, in which the bcc phase nu- 
cleates subsequently (see Fig. [T]). We have used the 
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FIG. 2: (color online) Nucleation map for PFC models: (a) 
IM-PFC; (b) 2M-PFC; (c) dependence on A at V^o = -0.25. 
The phase content obtained after 10^ time steps is shown: 
open triangle - liquid; square - (amorphous + liquid); circle - 
(amorphous + bcc); diamond - bcc; full triangle - amorphous. 
The heavy grey line stands for the stability limit of the liquid. 
Parts of the phase diagram are also shown. 



and ge order parameter to characterize the local structure 
[2|. With increasing undercooling, the nucleation rate 
of the amorphous clusters increases, leading to spatially 
nearly homogeneous transition at high undercoolings. In 
contrast, we have not detected any phase transition for 
more than 10^ time steps at e = —0.1598. These find- 
ings strongly indicate that crystal nucleation is enhanced 
by the amorphous precursor, and that direct bcc crystal 
nucleation from the liquid requires orders of magnitude 
longer time than via the precursor. This behavior ap- 
pears analogous to the role the non-crystalline precursor 
plays in colloids [TOl [TJ and simple liquids [8[ E] . 

We were unable to nucleate crystalline phases other 
than bcc in the IM-PFC model [Fig. [^Ja)]. Even in the 
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FIG. 3: (color online) Properties of the amorphous phase in 
the 2M-PFC limit from the ELE at e = -0.1: (a) Driving 
force for the liquid to amorphous transition (AcJamo-L). For 
comparison the driving force for fee freezing (Acj/cc-l) is 
also shown. The vertical lines indicate liquid densities at the 
(S) fcc-liquid and (MS) amorphous-liquid coexistences, (b) 
Nucleation barrier for the amorphous and fee phases at ipo = 
—0.21775. (a/cc - lattice constant of the fee structure.) Fits of 
W — Ar^-\-Br^ (solid lines) are compared with estimates from 
the classical droplet model using the equilibrium interface free 
energy and the driving force (dashed lines). 



stability domain of the hep and fee phases, the amor- 
phous phase formed in the time window of the simula- 
tions. Remarkably, this stayed so even in the 2M-PFC 
limit [Fig. [2|b)]. Interestingly, the amorphous phase ap- 
pears to coexist with the liquid, indicating a first-order 
phase transition between these phase, in agreement with 
the observed nucleation of the amorphous state. This 
suggests significant differences (e.g., in density) between 
the liquid and amorphous phases. Varying A at e = —0.1, 
we see a gradual transition from the IM-PFC behavior 
(liquid amorphous bcc) to the behavior seen on the 
2M-PFC side (liquid ^ amorphous) [Fig. ^c)]. (Com- 
parable results were obtained for constant cooling rates.) 

To investigate the lack of fee crystallization in the 2M- 
PFC model specifically designed to crystallize to the fee 
phase, we used the ELE at e = —0.1 for determining the 
free energy 7 of the fcc-liquid and amorphous-liquid inter- 
faces (Table I), and the driving force Acj for fee freezing 
and amorphization [Fig. [sj^a)] (see the methodology in 
fl5j). whose interplay determines the nucleation rate ^24j . 



In the density range of interest the fee phase is preferred 
thermodynamically [Fig. [sja)]. Between the fcc-liquid 
and amorphous-liquid coexistences, there is no driving 
force for amorphization, so fee freezing should take place. 
To evaluate 7, we have created equilibrium sandwiches 
(fcc-liquid-fcc and amorphous-liquid-amorphous), solved 
the ELE, and determined the grand potential emerging 
from the two interfaces. Comparable density changes 
were found at the amorphous-liquid and fcc-liquid tran- 
sitions (Table I), a finding consistent with the first order 
amorphization transition implied earlier. Remarkably, 
7am-L ^ ^'^^^y fcc-L (Table I). The nucleation barriers 
calculated for ipQ = —0.21775 using the classical droplet 
model, W{r) = (47r/3)r^Acj + 47rr^7, are shown in [Fig. 
[sjb)]. At this supersaturation fee nucleation is clearly 
preferable. Yet, we have not seen freezing even after 
1.9 X 10^ time steps. The liquid density beyond which 
Wfcc > Wam is ipo ^ —0.21. This is consistent with 
the finding that in the density range, where solidifica- 
tion could be observed at all (?/^o ^ —0.1962 for 10^ time 
steps), the amorphous phase nucleated. It appears that 
we cannot observe fee nucleation because of a technical 
difficulty: the time accessible for simulations is too short. 
Note that the dynamic EOM studies and the equilibrium 
ELE results consistently indicate separate time scales for 
changes of density and structure. 

For the PFC models the free energy barrier has a rough 
surface with many local minima that can be mapped out 
directly by using the ELE [15^. This provides an inde- 
pendent test of the previous computations. The results 
shown in [Fig. [sj^c)] indicate a reasonable agreement be- 
tween the predicted (dashed lines) and the directly eval- 
uated barriers (symbols + solid lines). This is reflected 
in the similarity of the interfacial properties evaluated in 
equilibrium and from fitting the droplet model (see Ta- 
ble I) . The minor discrepancy presumably originates from 
the fact that for such nanoclusters the classical droplet 
model is probably not very accurate. 

To rationalize the dominance of amorphous solidifica- 
tion in a substantial part of the phase diagram, we have 
evaluated effective pair potentials for the IM- and 2M- 
PFC models from the pair correlation function of the 
respective amorphous phases [Fig. [4]^a)] using Schom- 
mers' iterative method ^25^ that works reasonably for sin- 
gle component systems The potentials obtained are 
similar for short distances, and have a peak at ~ roA/2, 



TABLE I: Equilibrium densities ('0^'^) and interface free en- 
ergies (7) in equilibrium (eq) and from fitting to W{r) (fit) 
for the 2M-PFC model at e = -0.1. 
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eq 


jit 
Tx-L 


am 


-0.21885 


-0.21404 


1.79 X 10"^ 


2.34 X 10"^ 


fee 


-0.22139 


-0.21629 


2.76 X 10"^ 


2.80 X 10"^ 



Subscripts S and L stand for solid and liquid, respectively. 
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FIG. 4: (color online) Amorphous structure and effective pair 
potential in the IM- and 2M-PFC models: (a) g{r) of the 
amorphous structures, (b) Effective pair potentials evaluated 
from g(r) by Schommers' method. 



where ro is the radius at the main minimum of the poten- 
tial [Fig. |4jb)]. Remarkably, such potentials have been 
designed to realize monatomic glassformers, as a peak at 
roA/2 suppresses the close packed crystal structures [27 . 
Hence, we associate the evident difficulty to produce the 
fee and hep phases with this feature of the PEG effective 
interaction potentials. Furthermore, in the presence of 
multiple minima of the interaction potential coexistence 
of disordered phases is expected [2F, as indeed seen here. 

Summarizing, the PEG models display MS amorphous- 
liquid coexistence and first-order amorphization. In the 
cases accessible for dynamic simulations, the nucleation 
of the amorphous phase is faster than crystal nucleation. 
This leads to a separation of time scales for density and 
structural changes, as seen in other systems [9 . How- 
ever, some details differ: Such coexistence is unknown in 
the HS system, while the fee and hep structures are sup- 
pressed here. It is also unclear whether along the reaction 
coordinate specified in Ref. 8, the free energy landscape 
of the PEG models is similar to that of the LJ system. 

Gombining the results obtained for various potentials 
(LJ, HS, the present PEG potentials, etc. [8, 9, 27, 28 ), it 
appears that a repulsive core suffices for the appearance 
of a disordered precursor, whereas the peak at ~ 
correlates with the observed suppression of fee and hep 
structures, while the coexistence of the liquid and amor- 
phous phases seen here can be associated with multiple 
minima of the interaction potential. 
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